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Abstract 

o ; 

, A countable set of asymptotic space - localized solutions is constructed by the complex germ method 

in the adiabatic approximation for 3D Hartree type equations with a quadratic potential. The asymptotic 
parameter is 1/T, where T ^> 1 is the adiabatic evolution time. 

A generalization of the Berry phase of the linear Schrodinger equation is formulated for the Hartree 
type equation. For the solutions constructed, the Berry phases are found in explicit form. 

> ' 

F-» ; 1 Introduction 

O ■ 

The method of semiclassical asymptotics developed in [1-4] is used to construct and study the solutions 
of a many-dimensional non-stationary Hartree type equation with a nonlocal nonlinearity and the coefficients 
adiabatically varying in time. By adiabatic evolution of parameters we mean their relative small changes 
i in some characteristic time (adiabatic evolution time). In particular, for parameters periodic in time, the 
characteristic time is the period. 

In the time of adiabatic evolution, a quantum system returns to its initial state and the wave function 
i-Q gains only a phase factor. M. Berry has revealed [5] that the total phase contains, along with the dynamic 
d ■ part known from the Born-Fock adiabatic theorem [6-8] , a geometric additive component (geometric phase 
(GP)). This summand in the total phase is known as the adiabatic phase or the Berry phase. The adiabatic 
• • ■ phase is closely connected with the Floquet problem for systems of differential equations with periodic 
coefficients [9]. In quantum mechanics, geometric phases are well investigated for the linear Schrodinger 
equation [6,10-12]. In classical mechanics, a Hannay angle is introduced for nearly integrable Hamiltonian 
5h ' systems with adiabatically varying parameters. The Hannay angle is related to the Berry phase [13] if 
. i the Hamiltonian system corresponds to the considered quantum system. The Hannay angle is defined as 
an additional term to the dynamic contribution in the "angle" variable when the Hamiltonian system is 
described in terms of "act ion- angle" variables [14,15]. 

Geometric phases are observable [16-21] and they show up in various physical phenomena [22-24]. It 
is considered that in quantum calculations [25,26], intensely developing at present, the Berry phase might 
be used in some types of quantum gate, the so-called geometric gates (GGs) [27]. The last ones offer some 
advantages over the conventional (non-geometric) phase gates as the GGs have the greater fault tolerance. 
An example of an engineering embodiment of geometric gates based on the nuclear magnetic resonance is 
given in [28] . To build hardware for a quantum computer systems of cooled ions can be used [26] . In such 
a system, each ion caries a qubit, and the execution of a logic operation (gate) is governed by an external 
electromagnetic field created by laser radiation, a magnetic field source, etc. It should be noted that the 
Paul traps can be described by the potential of a harmonic oscillator [20,26]. One should also take into 
account the collective coupling of ions as the ion states depend not only on the external field, but also on the 
collective behavior of the ions [26] . A consideration of the coupling between parts of a system naturally leads 
to nonlinear models. An example of such a system might be the Bose-Einstein condensate (BEC) [29,30] 
models of which are used the Gross-Pitavskii equation [31,32]. This is an argument in favor of to study such 
models taking in mind their applications. 

The geometric phases in nonlinear systems were studied in [33,34] where the Berry phase was found for 
states describing a two-component BEC in an external field. In such systems, the Berry phase arises when 
the parameters responsible for the coupling of the condensate subsystems change slowly. 
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Here we continue our investigations started in [35, 36] where geometric phases were studied for a one 
dimensional Hartree type equation. The aim of this work is to find the Berry phase for a many-dimensional 
Hartree type equation of the following form: 



-iHd t + H„{R{t),*{t))\* = 



H„{R{t)^(t)) = H(R(t)) + xV(R(t), *(t)), 



H(R(t)) 



1 (u e r \ 2 p(t)((x,d + &g» k(t)x 2 
— \p--A{x,t)) + 



V(R(t),*(t)) = | / dy a(t)x 2 



2 2 

2b(t)(x,y)+c(t)y 2 ]\*(y,t)\ 2 , 



(1.1) 
(1.2) 



(1.3) 



A(x,t) = -[H(t),x\. 

Here, a(t), b(t), and c(t) are parameters of the nonlocal potential; k is a nonlinearity parameter; mit), kit), 
and p(t) are time-dependent parameters of the system, H (t) is an external magnetic force; (a, b) and [a, b] 
are, respectively, the scalar and the vector product in M 3 . The nonlinear Hamiltonian Ji^ depends on time 
via the set of parameters R(t)=(m(t), k(t), p(t), H(t), a(t), b(t), c(t)). 

The states with nontrivial geometric phases exist for quantum systems with Hamiltonians quadratic in 
variables and derivatives. Such Hamiltonians, in particular, include the external field of a harmonic oscillator 
and a uniform magnetic field [37]. 

The magnetic field reduces the symmetry of a quantum system, removing its degeneration, which in turn 
can result, for instance, in a normal Zeeman effect [38], etc. 

The study of geometric phases for nonlinear equations is a complicated mathematical problem since 
to construct the phase requires methods for solving many-dimensional nonlinear equations with variable 
coefficients in an appropriate class of functions. The well-known Inverse Scattering Transform method [39-41] 
is applied mainly to (1 + 1)D and (1 + 2)D nonlinear equations with constant coefficients, and only soliton 
solutions can be constructed in explicit form. The symmetry analysis [42-47] allows one to study systems 
possessing high symmetry, but the evaluation of a symmetry is inconvenient if the equation contains nonlocal 
terms. The nonlinear problems of the above class could be effectively solved by the method developed in [1-4] 
where asymptotic solutions (in some cases, exact solutions) are constructed for a nonlinear Hartree type 
equation which is the Gross-Pitaevskii equation with a nonlocal nonlinearity. 

This paper is organized as follows. In the first section, necessary designations and definitions are intro- 
duced for the Berry phase. In Section 2, a method of seeking an exact solution for a Hartree type equation 
with a quadratic Hamiltonian is briefly described. Section 3 presents a solution of the spectral problem for 
an instantaneous Hartree type operator. The solution is used to extract the dynamic phase from the overall 
phase. In Section 4, solutions of the Hartree type equation are constructed in the adiabatic approximation 
and the corresponding Berry phases are found. In Conclusion, the results and related problems are discussed. 



2 The Berry phase for a nonlinear equation 

In this work, we find the Berry phase using the approach developed in [37,48] for linear quantum equations. 
This approach is based on seeking an exact (or approximate) solution of the Cauchy problem for equation 
(01 : 

*U = Vv(x,i?(0)). (2.1) 

which then is expanded in an adiabatic parameter. Here the functions i/j v (x,R(0)) are determined by the 
spectral problem for the instantaneous Hartree type operator 

H„{R(t))MZ> R (t)) = E v {R{t))i> v {x, R(t)). (2.2) 

Assume that the parameters R(t) are T-periodic and slowly vary with time 4 . Following the adiabatic 
theorem in (linear) quantum mechanics [8] , we seek a solution of the Cauchy problem 1|1.1|) , l|2.1() in the form 

4 A system is assumed to be adiabatic if the following condition is fulfilled: 

■ T 

max — «1, (2.3) 

i— l,n R% 

where Ri are parameters of the Hamiltonian (see, e.g., [49]). 



t) = exp[i^(t)]^(f , Ji(t)) + Ol i ) . (2.4) 



T 



Then for t = T, taking into account that R(T) — i?(0), we have 



*(f,T) = exp[#,(T)]*(x,0) )• (2-5) 

Let us rewrite the phase 4> V (T) as 

I/ (T)=^(r)+7„(T) J (2.6) 
where 8 V (T) is the dynamic phase, which is fined by the relation 

T 

5 V {T) = ~J E u (R(t))dt. (2.7) 
o 

Following [5] (see also [10,12,36]), we call the phase 7i/(T) the adiabatic Berry phase or geometric phase for 
equation l|l.l|l . The Berry phase for the linear Schrodinger equation (x = in 1)1 .2)1 ) is determined by 

T 

lv {T) = i J (Mx, R(t))\i> v (S, R(t)))dt. (2.8) 
o 

Formula (|2.8|l is equivalent to (|2.6() in the linear case and requires an additional substantiation for nonlinear 
equations. 



3 Method of semiclassically concentrated functions 

Consider exact solutions for equation l|l.l|l following [1-4] . From these solutions we find solutions of the 
form (|2.4|l and obtain the Berry phase. 

Define the mean value for a linear operator A in a state ^(t) as 



(A(t)) =A*(t,h) = -__(*(t)|A|*(f)) = J dxV*(x,t,h)A(t)*(x,t,h), (3.1) 

where H*^)!! 2 = For the solutions *(i) of equation (fTTTf) we have 

d(A(t)) /dA(t)\ i 



dt 



(m) + iW~(tMt)),A(t)]-), (3-2) 



where [A, B]_ = AB — B A is the commutator of linear operators A and B. 

As in the linear case, we call equation 1)3. 2|l an Ehrenfest equation. From 1)3.2)1 with A = 1 it follows, in 
particular, that the norm of a solution of the equation 1)1.1)1 does not depend on time, i.e. 

||vl/(f,t)|| 2 = ||vl/(f,0)|| 2 = ||vl/|| 2 . 

Then it is convenient to go from the constant xtoa constant k = x\\^\\ 2 . 

Let us write down the Ehrenfest equation for mean values of the operator q: 



z 3 ,-(Az k Azi + AziAz k ); j,k, I = 1,6) ; z m = p m , z m+3 = x m , m = 1, 3; (3.3) 



where Az k = z~k — Zk9(t,K), Zk^i(t,h) — (zk). As a consequence, we obtain for the first- and second-order 
moments 

f = Jf) z (gy,R(t)) 7 

(3.4) 

I A 2 * = JSj zz (R (t )A 2 * - A 2 ^ z JR(t) J, = A 2 *. 



2m(t) 



E, f) px {R{t)) 



Sj pp (R(t)) Sj px (R(t)) 
%l x (R(t)) Sj xx (R(t)) 

/ 2m(t)p(t)c 
e 

-H 3 (t) 



2m(t)c 



H 3 (t) 
2m(t)p(t)c 

e 



Hi(t) 
2m(t)p(t)c 



(3.5) 



Sj xx (R(t)) = 



Am(t)c 2 



( -> 0/ s 0/ s 4m(t)k(t)e 2 
' H 2 {t)-H 2 (t) + 1 ' K ' 



V 



-H 2 {t)H 1 {t) 
-H 3 (t)Hi(t) 



H 2 (t)-H 2 (t) 



-H 1 {t)H 2 (t) 

4m(t)k(t)e 2 



-H 2 {t)H x {t) 



H 2 (t)-H 2 (t) 



-Hi{t)H 3 {t) 

H 2 (t)H 3 (t) 

Am(t)k{t)e 2 



J 



Here E is an identity matrix, B 1 is the transpose to the matrix B; ki(t) = k(t) + k(a(t) + b(t)), 
k(t) = k{t) + ka{i). 

The matrix of the second centered moment A2* has the form 



A 2 *(t, h) 



where 



^xx(t^) = ht kXl {t,K)hx3 = \\(Ax k Axi 

= h XkPl (t,h)\\ 3x3 



( cr* p (t,ft) (Tp X (ti h ) 
V a xp(^ h ) °"xx(<> S ) 

(Ap fc Aft))|| 3x3 , 



(3.6) 



3x3, 



k , / — 1 ; 3 j 



-(Ax k Api + ApiAx k ) 



3x3 



We call the system of equations (|3.4[1 the second order Hamilton- Ehrenf est system (HES) corresponding 
to equation Ijl.ip . 

The matrix of the second centered moments 5 A2 can be rewritten 



A 2 (t)=A(t)A 2 (0)A+(t), 
where A{t) is the fundamental matrix of the system in variations 

A = JSi„(R(t))A, A(0)=I. 



(3.7) 



(3.8) 



We denote by $(t, €.) = (Z(t, €.), A 2 (t, £)) the general solution of the system (|3.4|) . where £ are integration 
constants. 

Let us seek a solution of equation in terms of the anzats 



V(x,t,h) = <p(-^,t,Vh)exp j^(s(t,h) + (P(t,h), Ax)) 



(3.9) 



Here the function ip(^,t,\/h) belongs to the Schwartz space § in the variable £ = Axj\[h and regularly 
depends on \fh\ Ax — x — X{t, K). The real function S(t, ti) and vector function Z(t, K) — (P(t, H),X(t, H)) 
are to be determined. 



5 The subscript * can be omitted 



Let us expand the operators entering equation as Taylor series in Ax — x — x^(t,h), Ay = y — 

x^(t,h), and Ap = p — p^(t,H). Then equation Ijl.lfl takes the form 

{-ihd t +Sj(V,t) + (Sj z (V,t),Az) + ^{Az,S) zz (t)Az)}V = 0, (3.10) 



e 



2 



8m(t)c 



H 2 (t)xi(t,h) - ((H(t),x*(t,h))) t 

fe=l 

Sj z (*,t) = fi z fa(t,h),R(t)), % zz (t)=F> zz (R(t)). 



3 



Here ko(t) — k(t) + k{a{t) + 2b(t) + c(t)); the vector $j z (Q\z, R{t)) and matrix S) zz (R(t)) are given in i|3.5[l . 

Let us relate the nonlinear equation H3.10J1 with the linear equation that is obtained from H3.10|l by 
formal replacement of the mean values of the operators of coordinates, momenta, and centered moments of 
the second order g*(i, K) by the general solution g(t, <£) of the Hamilton-Ehrenfest system (|3.4|1 : 

{-ihd t + Sy(t, C) + (Sj z {t, £),Az) + ^(Az,f) zz (t)Az)}<S> = 0, (3.11) 



B*(t)-»a(t,c) 



We call equation l|3.11ll the associated linear Schrddinger equation. We can immediately verify that the 
function 

detCft) ) ex p{K^' e:) + ( ^' £) ' Aa?> + ^ (Aa? ' Q(<)Aa?> )}' (3 - 12) 

where 

t 

S(t, C) = / ((P(t, <£), £)) - C))«ft. (3.13) 



is a solution of equation (|3.11(l . Here Q(t) = B(t)C (t), and by B(£) and C(t) we have designated the "mo- 
mentum" and "coordinate" parts of the matrix solution of the system in variations (see [50] ) corresponding 
to the linear equation (|3.11|) : 

B = —Tt X p(t)B — Tt xx (t)C, (3 14) 

C = H pp (t)B + H px (t)C. 
Let us write down a matrix 

aw = ( ™ ) . 

The matrix A can be represented as 

A = (ax, 02,03), (3.16) 
where {a k } is a set of linearly independent vectors — solutions of the equation 

a h = Jf) zz {t)a k , /c=T73. (3.17) 

An operator 6 

a{t) = N a (a(t), J T Az) (3.18) 

is a symmetry operator for equation (|3 . 1 1|> if the vector a(t) is a solution of the system in variations (|3.17(l [50] . 
Let a(t) and b(t) be symmetry operators corresponding to the two solutions of the system in variations, a(t) 
and b(t), respectively. Then it is easy to verify that 

[a(t), &(<)]_ = ihN a N b {a(t),b(t)} = ihN a N b {a(0),b(0)}. (3.19) 

The last equality in Ij3.19|l is due to the Hamiltonian form of system H3.17|l . By braces we designate the 
skew-scalar product of the two vectors {a, 6} = (a, J J b). 

Assume that the system in variations (|3.17|l admits a set of three, linearly independent complex solutions 
a k (t) = (W k (t), Z k (t)) T satisfying the skew-orthogonal condition 

{a k (t),ai(t)}=0, {a k {t) 1 a* l {t)} = 2i5 kh k,l = T^. (3.20) 



6 we omit the dependence on h and £ 



Recall that the 6-vectors afe(t) and at(t), k — 1,3, serve as a symplectic basis in C„, and the 3-dimensional 
plane r 3 (Z(t, £)) with the basis afc(t) constitutes a complex germ on z — Z(t, £) [51,52]. 

Setting Nj = (ft) -1 / 2 in formula (|3.18|) . we compare the vectors a,* it) with the "creation" operators 
and the vectors Oj(t) with the "annihilation" operators a,j(t). Then, taking into account (|3.19|) . we obtain 
for the operators (t) , a,j (t) the bosonic commutation relations 

[a j (t),a k (t)}_ = [aj(t),a+(t)]_=0, [a^t), a+(t)]_ = 6 jk , j,k = T^. (3.21) 

Statement 1 Let cij(t) = (Wj(t), Zj(t)) be solutions of the problem (I3.17|) . I|3.20|l . Then the function 
<f>o(x,t,£) = |0, t, £) (|3.12|) is a "vacuum" trajectory- coherent state 

a,(t,£)|0,t,£) = 0, j = M, (3.22) 

/or i/ie associated linear Schrddinger equation (|3.11|l . 

Proof. Applying the "annihilation" operator Oj(i) to the function |0,i), we find that 

dj%t) = \0,t)[(Z j (t),Q{t)Ax) - (W 3 (t),Ax)]. 
It follows immediately that (|3.22l) is true as 

Q(t)Zj(t) = B(i)C'- 1 (t)£j(f) = Wj(t) 

according to the definition and properties of the matrix Q(t). 

Let us define now a countable set of states \v, t, <£) (exact solution of equation (|3 . llfi ) as a result of the 
action of the "creation" operators on the "vacuum" state |0,i, €.) H3.12(l : 

* v {2,t,£) = M,£) = -±=(a+(t,£)y\0,t,<£) = 
3 1 

= JJ-=(a+(t,C))"*|0,t,C), v = {v x ,v 2 ,vi), (3.23) 



fc=i v 

Notice that, in view of the explicit form of the "vacuum" state (|3.12|) and symmetry operators (|3.18|l for the 
functions <&„(:?, t) (|3.23|) constituting the Fock basis, the following representation is true: 



= N„$o(x,t)ne„(((R,t),t). (3.24) 

Here, Re„(((R, t), t) are many-dimensional Hermite polynomials, which are set by the matrix W(t) (see [53] 
for details), 

Re u (C(R,t),t) = (-lp^-2W(t)cJ -1, N v = " -L; (3.25) 

C(R,t) = ^(C^t^Ax, W(t) = C+(t)(C-\t)y. (3.26) 

The functions $„(x,i, £.) are solutions of equation l|3.10|l with a proper choice of £, such that the solutions 
$(t, €.) of the Hamilton-Ehrenfest system l|3.4|l coincide with the corresponding mean values 0$„(i, H, £) for 
the states (|3~2lt (see [54]). 

Let us designate this set of parameters by (t u ; then we have 

^ v (x,t,K) = $ v (x,t,~l I/ ). (3.27) 

The subscript v in shows that for each function ~^> v (x,t, ft) there exists its own set of the parameters <L U . 

The construction of solutions l|3.27|l uses solutions of two auxiliary ordinary differential systems: the 
Hamilton-Ehrenfest system <|3.4|l and system in variations l|3.17|l . 

For an arbitrary time dependence of the coefficients R{t), the solutions of these systems are unknown. 
However, if the system parameters depend on time adiabatically, we can seek a solution of the Hamilton- 
Ehrenfest system and system in variations in the form of a power series in the adiabatic parameter for which 
l/T is taken. These solutions allow one to construct a leading term of the asymptotics in the parameter 
1/T of equation l|l.lfl . For solutions of this type, the Berry phase can be found in explicit form. 



4 The spectral problem 



Consider the spectral problem (|2.2I) to state the Cauchy problem (jl.lH , (|2.1|) and find the dynamic phase 

by 

The solution of the spectral problem can be obtained from the non-stationary Hartree type equation , 
where the operator , H >e (R(i), vE'(t)) is replaced by Ti^R, ^(t)) with R = const. The solutions of equation 
1)1.10 . which have the form 

*(f,i) = exp{ - ~E v {R)t}i) v {x,R) (4.1) 



provide a solution of the spectral problem \2.2\ . Here i/j v (x, R) and E U (R) are eigenfunctions and eigenvalues 
of the instantaneous Hamiltonian TL >C {R 1 ?/v(i?)), respectively. 

The spectral problem is related to the stationary solutions (g(t, <Z) = 0) of the Hamilton Ehrenfest 
system (13.411 written for the stationary nonlinear Hamiltonian Ti^R,^). The solutions of l|3.4[l determine a 
stationary point g(R, € s ) in an extended phase space. Here by £ s we have designated a subset of constants 
separated from the set € by the condition of stationarity of solutions of the Hamilton-Ehrenfest system 
(see [36]). From 

JSj z (Q,R) = 



it follows that P(t, £ s ) = 0, X(t, € s ) = 0. 

The solutions of the Hamilton-Ehrenfest system for the second-order moments are obtained from solutions 
of the system in variations according to (|3.7|l . 

The linearly independent solutions of equation l|3.17|l normalized by the skew-orthogonality condition 
(|3.20|) can be written 



a v (t) = e*""%(fl), o 3 (i) = e^MR); 



fv( R ) 



ZUl a 

' {e v - i(-l)"e e ) 



\ \/2mu a 



, h{R) = 



m(p - _ \ 

e 

n 



J 



(4.2) 
(4.3) 



Here we have used the notation 



e n = (cos if sin 9, sin if sin 0, cos t 



e 



(sin ip, — cosy, 0), (4.4) 



eg = (cos <p cos 9, sin if cos 8, — sin ( 



n n (R) = (-i)v^M +UJa ( R ). o 3 (i?) = jA_^ (45) 



where 



eH , s I e 2 H 2 k ~ , . 

7 ? = 1,2; u c (R) = — ; oJ a (R) = \-r^2+--p 2 - (4-6) 
2mc V 4m c m 

The unit vector e„ specifies the direction of the magnetic field, and the set of vectors {e v , eg, e n } constitutes 
an orthonormal basis in M 3 . We suppose that the friquencies fli,£l 2 , ^3 do not satisfy the resonance relation 
litti + + ^3^3 — 0, where ii, I2, 13 are integers. 

The stationary solution A2(i?, £ s ) can be shown to exist if the solutions of the system in variations can 
be presented in the form 14.2(1 (see [55,56]). 

Note that, according to (|4.2|) . the matrices B{t) and C(t) can be written 



B(t) = B(R)A(t) = GB Q (R)A(t), C{t) = C{R)A{t) = GC (R)A(t), (4.7) 



where 



B (R) 



C (R) = 



( \fm{ p - ioJg) _ \pm{p - iug) 
m(ip + uj a ) ^/m(ip + oj a ) 







m(p — ifis) 



(4.8) 







V m ^3 / 

A(t) = diag{exp(^i(i?)t),exp(^ 2 ( J R)i),exp(ir2 3 (i?)t)}, 



G = (e„,ee, e n ) 



sm <£> cos y> cos w cos tp sin t 
- cos ip sin cos sin ip sin ( 
— sin 9 cos 



(4.9) 

(4.10) 
(4.11) 



It is easy to verify that the matrix G is orthogonal. 
In view of the relations 



fc=i 



det C(t) 



■ exp 



Q(i?) = GQ (R)G T , Qo{R) = dia,g(imuj a , imuj a , imQ 3 ) , 
1/2 



iV ri (detC(0)) =(7rfi)" 



•3/4 



(4.12) 

(4.13) 
(4.14) 
(4.15) 



we obtain from l|3.12fl the vacuum solution of the associated linear Schrodinger equation 

$ o(^ «,<£.)= exp { - ( - - £ ca xltXfe (J?, C a )t - - £ Q k t + - (x, Qo(R)x) J } 



(7r7i) 3 / 4 



fc=i 



fe=i 



i ( 5k H \ 

exp I - ^i^^2 ca x k x k (R,£ s )t+ -^2ci k t)}(t>o(x, R, 



fc=l fe=l 
Here x = G T x. In view of (|3.24|l and l|4.7() we find 

M -, ^ • / ~ 3 



(4.16) 



(^=) ^exp{-l(| ^ or^, (i?, £ s )f + ^ 7*O fc (v k + i) ^ jfofr R, € S )H U , 



where 7J„(£(i?)) are Hermite polynomials, which are set by the matrix V^(i?). Here 



(4.17) 



e(i?) = ^(C*( J R))- 1 A^=-^ 



w„ 

v/2n; 



X- 











1 

1 
1 



W(R) = -C+(R)(Co 1 (R)V-- 

For the considered Hermite polynomials, the following relation results from 14.7J1 : 

3 

He, (C(R, t), t) = exp [ - i ^ 0*1**1 H v (£(#), R) . 



(4.18) 



(4.19) 



k=l 



Evidently, to solve the considered spectral problem, it suffices to know only the sub-matrix of variances 
of coordinates instead of finding the complete matrix of the second moments. 



For the matrix of coordinate variances calculated for the states $„(£, t, <t s ), the following formula is valid 
(see, e.g., [50]): 

vxx = -[CWD^C+iR) + C*(R)D~ 1 C T (R)], D' 1 = ||(2^ + l)S kj (4.20) 



3x3 



Using (|4.20|l . we find 



= (vi + l) — + (v2 + l) — + ("3 + l)-ir- ( 4 - 21 ) 

' » 2 / moja \ 2 / moo a \ 2 J ml h 



fc=i 

Then the solution of equation is 

\» 



*{x, t) = t, o = f i= ) _L eX p | - « J2(n k + n fc ) (v k + i^|i/„(e(i?))0o(x, fi, 

Therefore, the eigenfunctions of the Hartree operator (|1.2f) are 
and the corresponding eigenvalues are given by the expression 

3 - 

E u (R)=hY l (^k(R)+^k(R))(^k + ^)- (4-24) 



fe=i 



5 Adiabatic approximation and the Berry phase 

Assume that the evolution of a quantum system goes adiabatically. This implies that the Hamiltonian 
parameters slowly change with time [see <|2.3[l ] and, along with the "fast" time t entering the time derivative 
operator, a "slow" time s can be introduced which the Hamiltonian parameters (R(t)—R(s)) depend on. Let 
the "fast" and the "slow" time are related as 

s = t/T, (5.1) 

where T is the evolution period of the system. 

As mentioned above, to find a solution of equation (|l.lfl in the adiabatic approximation, it is necessary 
to solve the Hamilton-Ehrenfest system and system in variations accurate to the second order in 1/T. 

The Hamilton-Ehrenfest system can be written as 



-4, = JSj z (si*,R(s)), 

^A' 2S - Jfi zz (R(s))Aw - A 2 *Sj zz (R(s))J, Al = A 2 , 



(5.2) 



where a' = da/ds. We seek a solution of this system in the form 

Z(t) = ZW( S ) + ±2W(s) +o(±)> = A 2 0) (*) + ^ A 2 1} (*) + ( 5 - 3 ) 
and obtain 

X {Q) {s) = X«(s) = 0, PW(s) =pW(s) = 0. (5.4) 

As in the spectral problem, we obtain the solution of the Hamilton-Ehrenfest system for the second order 
moments using the solutions of the system in variations. 

Making the change of variables by (|5.1|l in the system in variations (|3.17|) . we obtain 



-a'(t) = JSj zz {s)a(t). 



(5.5) 



Let us seek a semiclassical asymptotic solution of the system (|5.5H in the form 

a k (t) = e « n *M++>WMt) + o(Je), (5.6) 

/fc(*) = /f (s) + ^ft 1) (s) ) k = T^. (5.7) 
Substituting (|5.6|) in 1)5. 5|) and equating the coefficients of equal powers of 1/T, we obtain 

(JSj zz (s)-i<f>' k (s))fl°Hs)=0, 

(JSj zz (s) -i^WJ/^W = fW + ^W/f'W- (5-8) 

Then 

*i(s)=n fc (a)=n fc (ii(s)), fc = l73 (5.9) 

f^(s) = f v (R(s)), ri 0) (s) = f 3 (R(s)), 77 = 1,2, (5.10) 

where the vectors f v (R(s)), fs(R(s)) and the functions f2fe(-R(s)) are determined by l|4.3|) and Q4.5J1 . respec- 
tively. 

Let us decompose the vectors fi (s) in the basis vectors fjP (s) and (s): 

/* 1) w = E^ m (^)/i 0) (^)+^ m ( S )/i o) *( S ). (5.ii) 



Then we obtain 



Note that 



or, in matrix form, 



^) = ^{fi 0) '(s)ji 0) *(s)}, (5.12) 
M _ {/1 0) ^),/1 0) >)} {/i 0) (^/r'(^)} 

afemlsj " 2(o fe ( s )-n ro ( s ))' 2 (n m ( s ) + n fe ( s ))- l5 - idJ 



otk m (s) = -a* mk (s), Pkm(s) = (3 m k(s) (5-14) 



A(s) = -A+(s), B(s) = B T (s), (5.15) 
where A(s) = ||afe m (s)|| and B(s) = ||/3fc m (s)||. By analogy with (|4.7|l . we write 

Bit) = (B( )( s ) + ifl«( s ))A(°)( s ,T) + O^), 

C(t) = (CW ( S ) + i(5« (*)) AW (*, T) + O (± \ (5.16) 



A(°)( S ,T) = diag{T$i(s) +0i(s),T$ 2 (s) + </. 2 (s),T$ 3 (s) + 3 (s)}. 



Here 



B(°)( S ))=S(i2(s)), C(°)( S ) = (7(i?( S )), 

B^s) = B(°)(s)A T (s) + B (0) *(s)B T (s), (5.17) 
C«(s) = C (0) (s)A T (s) + C (0) *(s)B T (s). 

The matrices B^ \s) and C^^s) are determined by the vectors l|5.10|l . Hereinafter we omit the argument 
s or t if this does not lead to confusion. 
Note that 

-(mp)' , l n , , n ,, -(mp)' , 1,., _ . ,, -(mp)' 



Q(t) = dia,g(imuj a (s),imu; a (s), imfi 3 (s)) + o(^\ ■ (5.18) 



For the matrices A and B we obtain 
/ an 




7 



B 



Here 



V 2(n 3 -ni) 







(e' n , (ie v - eg)) 



V 2(fl 3 + fii) 



(mp — imuj a y 

(e' n , (-ie v + eg)) 



2(Q 3 - fi 2 

(mp — imwa)' 
4muj 2 





2(ft 3 + VL 2 ) 



(e' n , (-ie v - eg)) 



Q 3 
2uj„ 



2n 3 



2uj„ 



2(fii-n 3 ) 



7 



2(fi 2 -n 3 ) 



(e' n ,(-ie v + e g )} \ 
(e' n , {-ie v - eg)} 



"33 



7 \ 

(e' n ,(-ie v + e 9 )) 



2(0 3 + 



7 



(4H e r e «)) 



2(ft 3 + f! 2 )™ 

(mp — i?7if2 3 )' 



4m^3 



20 3 



(5.19) 



the matrix elements an, a 2 2, and a 33 are imaginary, and they are determined from the next-order approx- 
imation. We do not give them in explicit form as these functions do not contribute to the leading term of 
the asymptotic expansion. 

To construct solutions of the nonlinear equation , we need the matrix of coordinate variances 



a xx (t) = ^ [C(t)D^C+(t) + CftD-Wit)] , 



and it suffices to know this matrix accurate to 0(l/T 2 

1 



^(t)= ( ri J(s) + -a^(s) + O 



1 

f2 



Here 



r (0) = » 

XX j 



fir-. 



(5.20) 



(5.21) 



(5.22) 



+&(° ) *{D- 1 A - AD- 1 }^ 1 + C^iD^B + BD- 1 }C ,(0)+ ] . (5.23) 

The leading term of the asymptotics depends only on the sum of the diagonal elements of this matrix 
3 



k=l 



h 1 h(mp)' 
muj a T 2m 2 io\ 



+ 



1 



l> 2 



1 h(mp)' 



1 



1 h(mp)' 



mfl 3 T2m 2 nl 



O 



1 

T 2 



rriLUa T 2m 2 u\ 

In view of this, a solution of the Cauchy problem p.ljl . (|2.1(l can be represented in the form 

1 



T 



where 



*^(f,t) =exp<j - -T I E v {r)dT + i lv {s) \^ v (x,R(s)) 



(5.24) 



(5.25) 



The function (|5.25|l is a solution of equation (|l.lfl in the adiabatic approximation. The quantities 



3 / 1 

E v (s) = hJ2(nj(s) + (u s + - 



7=1 



(5.26) 



J'=l 



3=1 



are eigcnfunctions of the instantaneous Hartree type Hamiltonian r H >c (R(s),ip^(R(s))), and 7^(s) is 

3 r / j^n - / \ 3 

2 / 

kc(t) 1 (m(T)p(T)Y , / 1 



i - 



2m(r)^( r ) 



3 = 1 

(m(r)p(r)y 



dT 



2m(r)w a (T) 



dr + [u 3 + - 



1 - 



hc(t) 
2m(r)0|(r) 



(m(r)p(r))' 



dr+ 



+ ^(^2-^1) y (e^r^e^r)) - (ee(r),e v (r))(iT. 



2m(r)fi 3 (r) 

(5.27) 



The functions fij(s), w a (s), and Ctj(s) are defined by (|4.5(l . 1(4.6(1 . and 14.22(1 . respectively. The evolution of 
the function 1)5.25(1 in the period yields 



-T I E v {s)ds + i 7u (T) \^(x,0). 



^ { °\x, T) = exp 

We now use ((2.6(1 and ((2.7(1 to define the dynamic phase 

i 3 

3 = 1 



(5.28) 



(5.29) 



and the Berry phase 

2 



7,C0 = £ 



i=i 



2mw„ 



1 



2w a V m 



dm + dp) + [V3 + 



1 



2mf2| ! 
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20,3 V TO 



-dm + dp 



(i>2 — VlJ 



c 



[HidH-2 — H 2 dH 1 ]. 



(5.30) 



Here, H±, H2, and Ho, are components of the magnetic field, H = \J H\ + fl| + i?| ; C is a contour in the 
parameter space. 

The Berry phase 1(5. 30|) of the nonlinear equation p. 1(1 differs from that of the linear Schrodinger equation 
(x = 0) by the presence of the parameter k instead of k and by an additional summand proportional to Sr. 
Let a = c in 1(1.3(1 . Then for (£l n i/Q ) 2 << 1 we obtain 



i=i 



1 

2co»o V m 



-dm + dp 



^3 



c 






1 




2^ 


&a\~ n 

1 x=0 





dm, + dp I + I v 2 — v\ 



Ho 



H(H* + H$) 



ftnl = \ , ^0 = 3 L_ n = \ p, uj Q = uj a \ = 



e 2 H 2 
~T2> 



Arrive m 



P 



(5.31) 
[H 1 dH 2 -H 2 dH 1 ], 

(5.32) 



The classical analog of the Berry phase is associated with the Hannay angles (see, e.g., [14]). Geomet- 
rically, the Hannay angles similar to the Berry phases. The Berry phase 7^ of a quantum system and the 
Hannay angles 6.; of the corresponding classical system [13] are related by 



®i = -h— — =-— , 1 = 1,3. 

al, dvi 



(5.33) 



Here Ii are quantized action variables and Vi are quantum numbers. The differentiation in 15.33(1 implies 
that Vi is a continuous parameter. In the nonlinear case, it is naturally associate the the Berry phase ((5.30(1 
with an analog of the Hannay angle by formula 1(5.33(1 : 



e? 



1 - 


KC 


J- 


— dm - 

\m 


HdpJ 


2mu> 2 


2uj a 


1 - 






— dm - 

^m 


^dpj 


2mto 2 


2^~ a 



Ho 



H(H( + H 2 
H(H 2 + H 2 



-[H x dH 2 - H 2 dH x ], (5.34) 
-\H\dH 2 — H 2 dHx], 



9f-/ 

When p = 0, the "Hannay angle" 0f is equal to the solid angle d(C) under which the curve C can be 
observed from the origin point in the three-dimensional parameter space (Hi, H2, H3). 

We have defined the Hannay angle in terms of quantum mechanics, since the nonlinear problem requires a 
special study of the "classical equations" corresponding to the nonlinear "quantum" Hartree type equation. 
In our consideration, the part of these classical equations is played by the Hamilton-Ehrenfest system .4|> , 
which has no Hamiltonian form with respect to the standard Poisson bracket. 

6 Conclusion 

A countable set of solutions has been constructed for the Hartree type equation ljl.2|l in the adiabatic 
approximation and the corresponding Berry phases have been found in explicit form. 

In the linear case (Jk = 0), the Berry phase is determined completely by the solution of the Hamilton 
system and the corresponding complex germ [37]. In the nonlinear case, the Hamilton-Ehrenfest system 
involves the Hamiltonian system and the complex germ. For the quadratic Hartree type operator l|1.2l) . the 
Hamilton-Ehrenfest system (|3.4|) breaks into the Hamiltonian system with a self-action (the first equation 
in ijlQJ) and the system in variations H3.8|l . 

It should be noted that, as in the linear case, the leading term of the asymptotics (|5.25|) in the parameter 
1/T remains an eigenfunction of the nonlinear Hamiltonian 1|1.2[) at every point in time and only gains a 
phase factor. However, the Berry phase for the nonlinear Hartree type equation cannot be calculated by 
formula (|2.8(l . This is because the linear principle of superposition was essentially used in deriving formula 
(|2.8J) . The Berry phase is obtained by selecting the dynamic part (by formulas l|2.6|l . from the total 

phase incursion in the period T. 

There exists a limit as Zc — > for every expression obtained, and the limits agree with results of linear 
quantum theory (see, e.g., [10,12,37]). 

The geometric Aharonov-Anandan phase [57] is a generalization of the Berry phase for the case of 
arbitrary cyclic states (*(T) = exp[$(T)]*(0)). 

The geometric Aharonov-Anandan phase for the linear Schrodingcr equation coincides with the Berry 
phase in the limiting case of adiabatic evolution (T — > 00) [12]. 

The generalization of the Aharonov-Anandan phase for cyclic solutions of nonlinear equations with 
unitary nonlinearity, to which class equation (|l.l(l also belongs, is given in [58]. 

In the case of adiabatic evolution, the limit of the Aharonov-Anandan phase apparently, depends on 
the type of the nonlinear equation and on the class of functions in which its solutions are sought. The 
relationship between the Aharonov-Anandan phase as T — > 00 and the Berry phase for solutions of the 
nonlinear Hartree type equation Ijl.lfl in the class of functions (13. 9[) can be a subject of further research. 

From the view point of practical applications in quantum computations, the non-Abelian Berry phase is 
important [59,60] which corresponds to a degenerated spectrum of the Hamiltonian. The method developed 
in the work can be easily generalized for non-Abelian phases. Note that, in contrast to the linear case, 
the occurrence of a non-Abelian phase is not related to the degeneration of the spectrum of the nonlinear 
operator, but it is related to degeneration of the spectrum of the corresponding associated linear operator. 
A complete investigation of this problem is beyond the scope of this work and requires special consideration. 
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